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THE SAS-D NUTATION CONTROL SYSTEM 
Javin M. Taylor 

NASA-ASEE Summer Faculty Fellow 


ABSTRACT 


A control law is developed for the SAS-D nutation control system. Nutation 

is removed in a sub-optimal manner with respect to fuel consumed, but attitude 
errors are minimized. 

The research performed consists of an investigation of nutation theory and 
nutation due to energy dissipation, nutation detection analysis, nutation control 
analysis, and error analysis. The resulting nutation control system uses an ac- 
celerometer or rate gyro to sense the nutation angle 0 which varies sinusoidally 
at the nutation frequency fi. The sensed nutation angle is compared with a 
threshold £ t . If the sensed nutation is greater, a thrust pulse of duration equal 
to the period of one spin cycle is initiated. Since the threshold 8 is set equal 
to the amount of nutation that can be removed by one thrust pulse of duration 
equal to one spin cycle, the spacecraft nutation is reduced to near zero. 


Error analysis indicates that the nutation sensed, the threshold, the amount 
of nutation that can be removed during one thrust pulse, and the time duration 
of one thrust pulse are all sensitive to the spacecraft spin frequency imparted 
by the Delta rocket. Thus, the nutation control system was designed to be a- 
daptive to the possible variations in imparted spin frequency. 
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THE SAS-D NUTATION CONTROL SYSTEM 


INTRODUCTION 

This report describes the Phase B study of the active nutation control for 
the SAS-D spacecraft. Following Grasshoff [ 1] , the theory of nutation control 
for a asymmetrical vehicle} such as the SAS-D, is developed. Then an engineer- 
ing analysis is performed which evaluates various nutation sensor schemes and 
nutation control schemes. Specifically, nutation sensing methods using an ac- 
celerometer, a rate gyro, or a pre-positioned rate gyro are developed. Then, a 
nutation control law is developed and an error analysis performed. 

Nutation control is required because of energy dissipation within a spin 
stabilized spacecraft spinning about an unfavorable inertia axis. The energy 
dissipation time constant is the critical parameter in designing the active nu- 
tation control system. Preliminary analysis by other investigators predicts an 
energy dissipation time constant for the SAS-D spacecraft of greater than 10,000 
seconds. However, since assuming a large time constant and experiencing a 
much smaller one can prove disastrous, this analysis assumes a time constant 
of 500-1000 seconds. This conservative design will still be valid if future 
analysis and simulation substantiate the large time constant. 


NUTATION THEORY 


The rate of nutation D and nutation angle & can be derived from the basic 
equations for a rotating body 

H = M (1) 


where H is the total angular momentum and M is the total applied torque. 
Furthermore, considering motion in a moving reference frame described by 


• • • 

®MOV ~ ^REL + M x ^MOV 

where « is the total angular velocity, and selecting the axis system to be the 
principle axes so that 


(2) 
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H = I 


= 1 w = 


I Cd 

x x 


I Cd 

y y 




(3) 


where I is the total inertia, the equation of motion becomes 

I w + U X (I w) = M. 


(4) 
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(5) 


Defining the spin axis of the SAS-D to be the x-axis and assuming s 0, 
the remaining differential equations can be written in the form 



Following Grasshoff [ 1] , and using complex variables notation, the differential 
equation for the transverse angular velocity can be expressed as 


Cd (t) + i Qtd(t) = N 


where 

/n Y /2 

“ (t > = a V <*> + (t). 


(8) 

( 9 ) 


2 



ft = ± (ft y ft a ) 1/2 (sign positive for SAS-D), 


( 10 ) 



( 11 ) 


Thus, if the solution for equation (8) is non-zero, nutation, consisting of a nuta- 
tion frequency ft and nutation angle 9 , exists. 

Following Timoshenko and Young [ 29 , the nutation motion can be described 
by a body cone moving around a space cone, as shown in Figure 1 • 



For the initial condition oo Q = (0), the torque-free solution of equation (8) is 

given by 


cu(t) = w 0 e 


-i ftt 


Substituting equation (9) into equation (12), 

ft., \l/2 


^(t) - % 0 + w *o e “ int 


( 12 ) 


(13) 
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(14) 


Assuming zero phase reference 


/ 0 \l/2 
(t) = A y sin f2 t + i A 


Ci) 


cos fi t 


where 


A ,= 


■ * *© 


n„ \ i 1/2 




zo 


(15) 




and letting ft y /Ci z ) 1/2 = k, then 


(16) 


A y = k A 2 


(17) 


From Figure 1 » for s m all angles the nutation angle 9 is defined to be 


6 


[HJ + H2]1/2 
= tan Q = —I 1 max 


H* 


which reduces to 


6 - 


(*>_ - n 


(18) 


(19) 


Returning to equation (8), the complete solution for the differential equation is 
given by 

= “o e' in<t ' t » ) + JL [l _ (20) 


Now, let the torque N be defined as 


N = N 0 e 


i A. 


( 21 ) 


where \ is measured ccw from the SAS-D y-axis* 
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For a torque applied for (t ; - t Q ) seconds, a> (tj) is given by 


- % e 


*i 0( t i “ t q ) % 


a\ 


i n 




( 22 ) 


and through trigonometric manipulation 


®(t 4 ) 


-iH( t “t 0 > 

ox e 




sin 


2 


i [\-fJ/2 (tft.)] 
e 


(23) 


The magnitude of &> (t i ) is given by 




-•Oft,-,,) 2 N 0 . rn(t,-t 0 )' 

/■> n 4- c ( n 

i C^-H/2 ( 

n 1 2 



(24) 


From equation (24), the magnitude of co (t 2 ) will be minimized with respect 
to k if the torque N 0 is applied 180° from eu e _1 ^ 1 1 " l o > . Thus , 


\ L_2-=7r-n(t l - t 0 ) 


(25) 


k '= n - 


“ to) 


(26) 


Furthermore, after obtaining the optimum torque angle from equation (26), the 
torque duration is maximized if 


i.e., 


- V 
2 


77 -/ 2 ; 


(27) 


(t, - t # ) = 1/2 T fl 


(28) 


where Tq is the period of the nutation frequency. Thus, maximum nutation is 
removed during each nutation cycle if the torque is aligned as in equation (26) and 
the torque duration from equation (28) is one-half a nutation cycle. 
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Returning to equation (24) let A&> be defined by 


% 


2N c 

n 


sin 



( 29 ) 


which for proper choice of ^ is the maximum change in the transverse angular 
velocity w(t). From Equation (19), the maximum change in the nutation angle is 


A <9 = 



(30) 


To su mmar ize, equation (29) shows the amount of transverse angular 
velocity that can be removed so that in equation (22), ^(tj) is less than o> 0 . 
Furthermore, since energy dissipation exists in the spin stabilized SAS-D, the 
transverse angular velocity tends to increase and must periodically be reduced 
through the use of nutation control. Increase in nutation due to energy dissipa- 
tion is discussed in the next section. 

Table 1 presents the analysis data for the SAS-D nutation control system. 
Lines 1-11 are the result of the material developed in this section. Nutation 
control will possibly be required during three parts of the spacecraft trajectory. 
These parts are: segment A, the transfer ellipse; segment B, after apogee 
burn; segment C, during synchronous orbit trim. Hence, the three columns in 
Table 1. 


NUTATION DUE TO ENERGY DISSIPATION 

Spin stabilization of the SAS-D spacecraft about an axis other than that of 
maximum inertia results in a possibly unstable stability condition. Energy 
dissipation due to fluid motion in heat pipe and fuel tanks causes the nutation 
angle to increase until the stable condition of spin about the, maximum inertia 
axis is reached. To counteract this condition, the nutation control system 
must add energy to compensate for the energy dissipated. 

It is generally accepted that for small perturbations the increase in nuta- 
tion angle due to energy dissipation is exponential and of the form 

6 (t) = e Q e t/r (30) 
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Table 1 

Analysis Data for SAS-D Nutation Control 


Item 

No. 

Parameter 

Segment A 

Segment B 

Segment C 

Comment 

1 


70.9 

62.4 

60.3 

Fig. 5.8-1, Phase A 

2 


211.8 

156.3 

151.0 

Report, units axe 

3 


216.7 

160.8 

155.5 

slug ft 2 

4. 


0.6884oj x 

0.6296^ 

0.6305 o> x 

Calculated with &> K 
as a parameter 

5 

6 


0.6502 
0.6690 co x 

0.5839 

0.6063o; x 

0.5833 w x 
0.606 4 co x 

: '.7: J 

k. 

1.0290 

1.0384 

1.0397 


8 

CO 

X 

1.0 rev/sec 

1.0 rev/sec: 

1.0 rev/sec 

Range is 1.0 rev/sec 
±20% 

9 

n<i.o) 

0.6690 

0.6063 

0.6064 

Units are rev/sec 

10 

T n (1.0) 

1.4948 sec 

1.6493 sec 

1.6491 sec 


11 

No 

0.0673 

0.0895 

0.0924 

Units are lb-sec 2 -ft. 

12 

Cr 

y 

0.3347 

0.3992 

0.3993 

V 1 * 71 * 

13 


0.3272 

0.3881 

0.3878 

O-, = I. /I, 

14 

r 

y 

31.552 sec 
8.76 hr 

12,359 sec 
3,44 hrs 

11,963 sec 
3.32 hrs 

Z X Z 

Energy dissipation 
time constant 



34,968 sec 

14,144 sec 

13,713 sec 

15 

T 

t 

9.71 hrs 

3.93 hrs 

3.81 hrs 


16 


30.4° 

19.1° 

19.1° 

1 

Nominal accelerom- 
eter position 

17 

l>J 

3.47 r # 

3.97 r # 

3.97 r # 

l 

• ^ , 

Accelerometer output 
magnitude 

18 

O 

PC 

-e- 

30.4° 

19.1° 

19.1° 

Nominal rate gyro 
position 

19 


2.94 # 

3.49# 

13.49# 

1 

Rate gyro output 
magnitude 

20 

^RG 

1 

! 

-59.55° 

-70.84° 

-70.84° 

Nominal phase for pre- 
positioned rate gyro 

21 


2.14 # 

2.57# 

2.57# 

Pre-positioned rate 
gyro output magnitude 
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Table 1— (continued) 


Item 

No. 

Parameter 

Segment A 

Segment B 

Segment C 

Comment 

22 

A 6 

Ct) 

X 

-0.7715° 

-0.7839° 

-0.8183° 

Thrust duration is one 
spin cycle 

23 

T 

* 

2303 sec 

1152 sec 

1152 sec 

Time between thrust 
pulses 

24 

H n 

1.833 lbs. 

0.1 lbs. 

0.76 lbs. 

Hydrazine reg T d for 
nutation control 

25 


-0.8752° 

-1.0490° 

-1.0813° 

Maximum 8 that can 
be removed 

26 

Max 6 

CD X 

515° 

237° 

247° 


27 

Max e opt 

585° 

317° 

327° 


28 

T 

1000 sec 

500 sec 

500 sec 

Conservative estimate 
of energy dissipation 
time constant 


where t is the time constant of energy dissipation in the spacecraft. 


In "Phase A Report for Small UV Astronomy Satellite," r is calculated to 
be on the order of 20,000 seconds* However, because of the inadequate 
theoretical and experimental knowledge at this time, a time constant of 20,000 
seconds is viewed with great skepticism. Thus, the energy dissipation time 
constant used in this report will be very conservative and on the order of 
500-1,000 seconds. 

Following Fedor T s memo [41* which uses the energy sink approach, the 
time constant due to the heat pipes for small nutation angles is given by 


where 


4 I T (l-o) 

T = 1 

m Cl) 2 Z? CT ^ 

3 1 

a - ratio of rotational to transverse interias for 
symmetric body. For SAS-D, cr = I A 

X T 

I T - transverse inertia for symmetric body. 
m ~ mass of fluid in heat pipes. 


(31) 
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, Z zf - summation of heat pipe locations. 

o> s - spacecraft spin. = &j x 

Using the parameters given in the Phase A report [3], page 5-69, for m and 
1 z? ,r ranges from about 12,000 to 32,000 seconds, (See items 12-15, 

Table 1.) 

Even if the energy dissipation due to fluid movement in the spacecraft is 
well understood, the time constant r varies greatly for small changes in the 
inertia ratio o - . Since o- will probably change during the course of the SAS-D 
program, an analysis of variation in r due to cr is warranted. Thus, the Taylor's 
series expansion of equation (31) about the nominal inertia ratio a Q , with higher 
order terms ignored, yields 



By allowing the approximation 

3- 2 cr 0 - 3 ( 1 - <t q ) (33) 

and simplifying, the final result is 



Thus, for example, a change in a from 0.330 to 0.397 (see Table 1, line 12, 
segments A and B), that is, 20%, results in a new time constant approximately 
40% of its original value. (See Table 1, line 14, segments A and B.) This fact 
reiterates the need to assume a conservative time constant at this point in the 
program. 

The increase in nutation angle described by equation (30) is valid for small 
values of 0 m Alfried [5l indicates that the exponential increase in nutation angle 
changes to the form 
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(35) 


5(t)=e o yT3: 

where r' is a different time constant from r. Here the nutation angle builds ud 
faster than in equation (30), p 

Furthermore, Alfriend gives the relation for the crossover from behavior 
described by equation (30) to equation (35) by 


b cr 0 


(36) 


where 

cr - ratio of rotational to transverse inertia 
Tj ~ damping coefficient 
b - heat pipe height to radius ratio 
H - nutation frequency. 

Without going into the detail, rough calculations indicate that 6 can build 
up about 15° before equation (35) is applicable. 


NUTATION DETECTION 


The SAS-D nutation can be detected using several methods. In this report 
nutation detection using an accelerometer or a rate gyro will be discussed. 

Accelerometer. To detect nutation using an accelerometer consider the 
acceleration at a point located P = col [x„ y 0 z ] from the spacecraft c.g 
The acceleration experienced at this point is 


a-wx p + tO X (tOX p) 


(37) 


which reduces to 


~ “ 


— — 


— 

£l 

* 


CO 2 . — CO V 

y 0 * ■'o 


" ("« y 0 - "y x o) + ", K Z o - ", V 

a 

y 


to X - CO z n 
* 0 X o 

+ 

", C"y x 0 - ", y 0 ) + ", (" y z o - ", y 0 ) 

- _ 


to x y n - cb x n 

* 0 y 0 


", X 0 - ", Z 0 ) + "y y 0 - " y Z„) 


( 38 ) 
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For an accelerometer mounted at x Q =0 with its input axis parallel to the 
principle x-axis, 

a * - ( W y - "«) y 0 + (V Z ^X + ^y) Z 0 < 39 ) 

Using the results of equations (9) and (14) that 

w y = ^ sin H t 
= A 2 cos H t 

equation (39) becomes 

a x = (Ay ^x + A 2 n)y 0 sin H t + (A e cn x + A y 0) z Q cos fi t (40) 

Again, following Grasshoff [1] and using phase notation, let 


a* = a Q Q s in (Q t + <f> ), 


(41) 


a 0 = C{ ( k ", + fi) y 0 > 2 + {(", tkfi) z 0 } 2 ] 1/2 (a, - n), (42) 


and 


4> 


(oj x + k H) z Q 
(k + n) y 0 


(43) 


Here <fi is measured ccw from the positive principle y-axis. Thus, from equa- 
tion (41) the acceleration sensed is proportional to the nutation angle £ and is 
sinusoidal at the nutation frequency rate, with phase shift 0. 


The phase in equation (43) depends upon the position of the accelerometer. 
Thus, it is possible to place the accelerometer so that if the thruster is actuated 
at some phase with respect to the maximum acceleration sensed, the nutation of 
equation (24) will decrease. Let equation (41) be modified to the form 


= a 0 & sin (fl t + /x) 


(44) 
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and 


ii = <t> + 0 + P 1 


(45) 


where 

0 - total phase shift due to electrical hardware. 

/3 ~ phase advancement to compensate for position of thrust centroid. 

Now, for optimum nutation removal, y. should be set equal to the optimum 
phase in equation (26). Hence, 


M \ - TT + 0 


(^) 


(46) 


From equation (45), 




<^+0+/5 t =X-7T + Q 

T \ 2 


(47) 


0 = A. — 77 — 0 + 


n 




(48) 


Substituting equation (48) into equation (43) the position of the accelerometer is 


k oj x + 0 


co + k n 

X 


t an <p 


(49) 


Thus, the angular position of the accelerometer is given by 


tan 0 A _ (50) 

y o 

where 4> k is measured ccw from the positive principle y-axis. 

For the SAS-D vehicle, nutation control will be accomplished by the two 
5-lb thrusters, with three foot moment arms, mounted on the plus and minus 
y-axis. Thus, k - ±tt/2. Furthermore, assume that 0 , fi - 0 in equation (48). 
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Also, assume that the applied thrust pulse will be the duration of one spin cycle, 
i*e., one second for a nominal « x =1.0 rev/sec. From equation (48), for 

\ ~ + 77 /2 


0 = (0 - 0 . 5 ) 77 


(51) 


where ft is in rev/sec. For X — - ~n /2, 


0= (n- 1.5) TT 


(52) 


Equations (51) and (52) show that the phase angles differ by 180° depending 
upon which thruster is fired. However, the phase angle of equation (51) can be 
used if the polarity of the accelerometer output is reversed for the alternate 
thruster. The nominal accelerometer location 0 A can be determined using 
equations (49) - (51). This angle is shown in line 16 of Table 1. Because of 
the value of k, 0 A - 0. 

Equation (42) gives the expression for the accelerometer gain. Since k (see 
line 7, Table 1) for the SAS~D spacecraft is close to one, equation (12) can be 
approximated by 



- ft) [yj + zj] 1/2 ‘ 


(53) 


Letting r = [y 2 + zj] 1/2 » the output of the accelerometer is given by 


\ = ~ ft 2 ) r 0 sin (fl t + 0), (54) 

where r is the radial placement of the accelerometer . This result is given in 
Hue 17, Table 1. 

Hate Gyro. To detect nutation using a rate gyro consider a rate gyro 
mounted at an angle 0 RQ from the principle y-axis. The measured angular 
rate is given by 


"ii = w y cos 0 RQ + sin 0 RG . 

Using the results of equations (9) and (14), equation (55) becomes 


( 55 ) 
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(56) 


~ \ sin 0 t cos c^ G + A z cos Q t sin <£ RQ . 

Furthermore, using the results of equations (17) and (19) let 

a M = % 9 Sin t + 0) ( 57 ) 

where 

% = [k 2 + 1] (co x - 0) (58) 

and 

tan ^Kc] ( 59 ) 

Now, by changing equation (57) to 


= 6 sin 1 + M) (60) 


the phase requirement analysis for the accelerometer is directly applicable to 
the rate gyro. Thus, equations (45) - (48) and (51) apply. Since k * 1, <t> - 0 

from equation (59). RG 

For the rate gyro, the gain from equation (58) is approximately /2(co - n). 
Thus, the output of the rate gyro is given by x 


= /2 (cu x - 0) e sin .(fl t + 0). (61) 


4 > = ta 


n-‘Ti 

_k 


The mounting angle and output are given in lines 18 and 19, Table 1. 

Pre- P ositioned Rate Gyro. Nutation can also be detected by a rate gyro 
positioned at any angle if the phase of the output is adjusted electronically. 
Considering that the y-axis rate gyro is used, equation (55) becomes 


CO 


M 


CO 

y 


(62) 


and equation (60) reduces to 
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£<j u = - k (oj x - O') 6 sin (H t - 77/2 + ij.). 


Since 0 = 0° for the y-axis rate gyro, from equation (45) 


or 


-7r/2+At=0+^ + y5 T 


fJ- ~ rr / 2 + 0 + / 5 t 


Setting equation (46) and (60) equal 


K -tt + Cl L_ — 2.^ = tt/ 2 + ^ + Aj. 


4 > = \ - 5/2 77 + n 


Assuming /? T = 0 and (t x - t Q ) = 1 second, for \ - + 77 / 2 , 


(63) 


(64) 


(65) 


( 66 ) 


'Aro = ( n - !>*■ 


for A = - 77 / 2 , 


(67) 


0r G = ^77. (68) 

Equations (67) and (68) differ by 180°. The phase angle of equation (67) can 
be used for both thrusters if the polarity of the rate gyro is reversed. The 
phase angle in equation (67) is given in line 20, Table 1. The output of the rate 
gyro, from equation (63) is given in line 21, Table 1, 


NUTATION CONTROL 

In discussing a nutation control system it is worthwhile to consider the 
effects that nutations will have upon the SAS-D mission. Upon entering the 
transfer orbit the Delta vehicle spin stabilizes the SAS-D spacecraft so that 
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the apogee burn can be oriented properly for synchronous orbit insertion. 

Thrust misalignment and nutation tend to divert thrust so that the Av imparted to 
the spacecraft for synchronous orbit insertion is not maximum. Considering 
the misalignment angle very small, the loss in Av due to a nutation angle Q is 
given by 


SAv^-!!av. (69) 

2 

Since synchronous orbit trim requires approximately one pound of hydrazine 
for each 3 m/sec correction in Av , minimizing the nutation angle 6 , in general, 
minimizes the hydrazine budget required for trim. 

In some spacecraft the nutation angle must remain small so that ejection of 
the apogee kick motor will not strike the spacecraft structure. This is not a 
consideration for the SAS-D since the apogee kick motor will not be ejected. 

The most significant reason for requiring a small nutation angle is apparent 
from equation (30). The smaller the initial nutation, the greater the time elapsed 
before the nutation builds up to a given level. Hence, nutation removal is re- 
quired less frequently and hydrazine is conserved. 

Grasshoff [ 1 ] describes a nutation control method in which the sinusoidal 
sensor output is compared against a threshold. The thruster is fired for the 
duration that the sensor output is greater than the threshold. Nutation is re- 
moved by a series of thrust pulses until the nutation is below the threshold. 
However, if a mi nimum thrust pulse length and/or secondary threshold are not 
employed, near the threshold the thrust pulse is so short that nutation removal 
and build up are equal and hydrazine is wasted. Furthermore, if the thrust 
pulse duration is not an integral multiple of the spin frequency period, the total 
momentum vector orientation is altered. 


For the SAS-D spacecraft a more straightforward nutation control can be 
employed. In this case, each time the thruster is fired, it fires for the period 
of one spin cycle. Equations (29) and (30), repeated here for convenience, 
determine the nutation angle removed by one thrust pulse. 


A co 


0 



sin 



A $ 


~ Aa) o 
"x- 0 


16 



Thais, for (t l - t 0 ) = 1 second, which is the nominal period for one spin cycle, 
A 6^ can be calculated. See line 22, Table 1, 

A A 0^ of about 0.8° can be removed by a single thrust pulse of one spin 
cycle duration. Hence, the nutation control can be designed with a threshold 
=0.8°. Whenever the sensed nutation angle exceeds this threshold a thrust 
pulse of one spin cycle duration will be fired. 

Once the type of nutation control described above is established two ques- 
tions must be answered: What is the hydrazine budget required, and is it near 
optimal? To answer the first question equation (30) can be re-written as 



Assume that the nutation threshold is set at # T = 0.8° and that the thrust pulse 
removes 90% of the nutation. Thus, & r = 0.8, =0.08 and their ratio is 10. 

Also consider that the transfer ellipse (segment A in Table 1) may last until 
3rd apogee, which is approximately 28 hours. Furthermore, assume a conserva- 
tive dissipation time constant r = 1000 seconds. From equation (70) the time 
between thrust pulses is T N = 2303 seconds (see line 23, Table 1). Thus, 44 
seconds of thrust are required for the 28 hour (maximum) transfer ellipse. As- 
sumping the low duty cycle specific impulse of hydrazine to be 120 Ibf-sec per 
lbm (see page 5-67, Phase A report [ 3} ), the weight of fuel required for a 5-lb 
thruster is 1.833 pounds (see line 24, Table 1). 

Lines 24 and 25 of Table 1 also show the time between thrust pulses and 
hydrazine required for segments B and C. Here it is assumed that r = 500, 
segment B lasts for one hour and segment C lasts for eight hours. 

To answer the question about the optimality of the proposed nutation control 
system consider equation (28) which shows that maximum nutation angle is re- 
moved if the thrust pulse duration is one-half the period of the nutation frequency. 
A comparison of the nutation period (line 10, Table 1) with the one second period 
of the spin frequency indicates that the thrust pulse is applied for more than half 
the nutation period. Thus, the nutation control system is sub-optimal in terms 
of hydrazine fuel consumed. 

The optimal nutation removal can be determined from equations (29) and 
(30) with (tj - t p ) = 1/2 Tq. These results are shown in line 25, Table 1. 
Comparing the proposed thrust duration with the optimum thrust duration 1/2 T 
(line 10, Table 1), it appeals that about a 34% fuel waste occurs. However, 
since a thrust duration of pne-spiii cycle causes no change in the momentum 
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vector, and since the sub-optimal hydrazine budget is small, the proposed 
system is recommended. 

Another advantage of the proposed nutation control system results from the 
behavior of the hydrazine thrusters. Considerable difference occurs between 
cold and hot hydrazine thrusters in terms of specific impulse and thrust pulse 
shape. With the proposed nutation control system, time between thrust pulses 
is at least 1152 seconds (see line 23, Table 1). Thus, the thrusters can be con- 
sidered cold at each firing and the nutation control system does not need to be 
adaptive to hot thrusters. 

Figure 2 shows a block diagram of the proposed nutation control system. 
The phase of the accelerometer or rate gyro output is set (by physical position 
or electronically) so that a thrust pulse of one spin cycle duration is activated 
if the nutation senses exceeds the threshold. The threshold is a critical param- 
eter and depends upon the amount of nutation which can be removed during one 
spin cycle. Thus, the spacecraft spin must be sensed to determine the thrust 
pulse duration and threshold. 


THRUSTER 



Figure 2. Nutation Control System. 

Figure 3 shows the trajectory of the spin axis oj during nutation buildup 
and removal. The spin axis position spirals out until the threshold 6 is 
reached. A properly phased thrust pulse is initiated which drives the spin axis 
to near zero. 

A final question concerning the nutation control system needs to be answered. 
Considering the energy dissipation model of equation (30) and (31), how great 
can the nutation angle become before the nutation buildup exceeds the amount of 
nutation that can be removed? Equation (30) can be rewritten in the form 

@ Q - A 6 - & 0 e t/T 
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Figure 3. 


Trajectory of Spin Axis During Nutation Buildup and Removal. 


where A# is the amount of nutation that can be removed by a thrust pulse and t 
is the nutation period. Then, the solution for & 0 = max 9 gives the nutation 
angle at which nutation due to dissipation equals nutation that can be removed. 
These results are given in lines 26 and 27, Table 1, for a thrust pulse of one 
spin cycle duration and for an optimum thrust pulse. Thus, for practical nuta- 
tion angle limits of 15° - 30°, nutation removal is considerably greater than 
nutation due to energy dissipation. 
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error analysis 


In this section of the report it is worthwhile to look at an error analysis of 
the nutation control system. The two most important variables to investigate, 
in terms of error, are the nutation build-up due to energy dissipation and the 
amount of nutation removed per thrust pulse. 

The nutation build-up due to energy dissipation will be analyzed first. The 
basic relationships are given in equations (30) and (70) and are repeated here 
for convenience. 


0 (t) = 6> 0 e t/T 



The time T N between thrust pulses depends upon the energy dissipation 
time constant r and the ratio of the threshold & r to the residual nutation 6 Q 
remaining after a thrust pulse. From equation (35), largest changes in r result 
from errors in the inertia ratio. Changes in T n due to changes in the amount 
of nutation removed per thrust pulse can be determined by a Taylor’s series 
expansion and retention of the linear terms. The final result is 


T = T fl-3 Scr 

" l L cr 


For 90%, or better nutation removal per thrust pulse, i n (6> T /e ) is greater 
than 2.3. Hence, a 10% error in the threshold setting results in about 5% de- 
crease in the time between thrust pulses. Thus, 5% more hydrazine is required. 

A larger effect occurs due to errors in the inertia ratio cr * As equation (71) 
shows, a 10% change in inertia ratio results in a 30% decrease in the time 
between thrust pulses. 





1 _ 

i $ e 



/6>\ 6 



ln (f ) 




\ o 



(71) 


Nutation removal will be discussed next. The basic relationships given in 
equations (29) and (30) are repeated here for convenience. 


A CO 


o 
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Changes in the amount of nutation removed can result from the thrust pulse 
duration (tj - t 0 ), the spin frequency a; x , and the thrust magnitude N 0 • Since 
thrust pulse duration is set as the period of one spin cycle, thrust pulse dura- 
tion errors can be analyzed as spin frequency errors. Thus, the principal 
error sources are and N 0 . 

Combining equations (29) and (30) and expressing Cl as Cl (co x )co x (see line 
6, Table 1); the expression for A 6 is given by 


- 2 N 0 sin 
A $ =. 

(n (*>„ 




o> ) CO ( t. - 

X' X ' 1 
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) - n 2 co 2 


(72) 


The change in A0 due to changes in <o x can now be determined from a Taylor T s 
series expansion and retention of the linear terms. Thus, 


A -0 i A ft, 


• t 

[i "(i _nK) (V^) cot 


(73) 


Letting co x (tj - t 0 ) = 277, the final result is 


A e 2T A e . 


1 - 1.6 


S co 


co. 


] 


(74) 


Additional error analysis with respect to N 0 yields 



(75) 


From equations (71) and (75) it is evident that the time between thrust pulses 
T n is sensitive to changes in the threshold and that the nutation removed per pulse 
A^ is sensitive to changes in the spin frequency co x . Since the nutation control 
law is based on the fact that the threshold <9 T = A# , the nutation control must 
be adaptive to the spin frequency. 
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An error analysis must also be performed upon the output of the nutation 
sensor (accelerometer or rate gyro) since the output depends upon the spin 
frequency. From equations (54) and (58) it is evident that the amplitude of the 
sensor output varies with changes in the spin frequency. For the accelerometer 
the output change varies as twice the spin frequency change. For the rate gyro 
the variation is direct. 

To summarize the results of the nutation control system error analysis, for 
the rate gyro sensor a 10% increase in spin frequency causes a 10% increase in 
the sensor output for a given nutation angle 6 . In addition this change also 
causes a 16% decrease in the proper level of the threshold 6> T . The result is a 
threshold that is too low for the nutation removed by the thrust pulse and over 
correction occurs. 


On the other hand, a 10% decrease in spin frequency causes a 10% decrease 
in the sensor output for a given nutation angle e . In addition, this decrease also 
causes a 16% increase in the proper level of the threshold 0 . The result is a 
threshold that is too high for the nutation removed by the thrust pulse and under 
correction occurs. 

In either case, the ratio (& t /& q ) in equation (70) is reduced and the time 
between thrust pulses T n is diminished. The result is additional fuel consump- 
tion. Thus, the nutation control system, as shown in Figure 2, must be adaptive 
to the spin frequency. > 

Once system performance degradation due to spin frequency error has been 
eliminated the remaining error source is the thrust magnitude N 0 . Considering 
the thrust magnitude to be known to within 10%, thrust errors will reduce the 
ratio (0 T /d 0 ) in equation (70) to 10 which is, in fact, the ratio that the nutation 
control system design is based upon. Thus, this error has been properly treated* 


CONCLUSIONS 

This report has described the analysis for the SAS-D nutation control sys- 
tem. The results are a control system design, predicated upon a conservative 
energy dissipation time constant, that provides properly phased thrust pulses 
of one spin cycle duration. Nutation is removed in a sub-optimal manner with 
respect to fuel consumed, but attitude errors are minimized. The analysis 
considers the possibility of sensing nutation via an accelerometer, a rate gyro, 
or a pre-positioned rate gyro. 


The results of this study should be refined when a better knowledge of the 

spacecraft energy dissipation and hydrazine thrust pulse characteristics become 
available. 
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